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Periodic  Orbits  about  an  Oblate  Planet 
by 

Richard  B.  Barrar 

ABSTRACT 


Periodic  solutions  of  the  equations  of  motion  of  a  satellite  of  an  oblate 
planet  are  investigated,  using  methods  developed  by  Poincar£  (1892).  The 
treatment  covers  both  critical  and  non-critical  angles. 

INTRODUCTION 

Much  work  has  recently  been  devoted  to  the  motion  of  a  satellite  of  an 
oblate  planet;  see  Brouwer  (1959),  Garfinkel  (1959),  Kozai  (1962).  However, 
very  little  attention  has  been  paid  to  the  existence  of  periodic  orbits. 
McMillan  (1910)  discussed  periodic  orbits  about  an  oblate  planet.  However,  all 
his  orbits  reduce  to  circular  orbits  in  the  unperturbed  case.  In  the  present 
paper,  we  apply  methods  developed  by  Poincare  (1892)  to  investigate  periodic 
orbits.  We  cover  both  critical  and  non-critical  angles.  It  is  found  that  the 
periods  are  completely  different  at  critical  and  non-critical  angles. 

Background 

In  order  to  make  the  presentation  as  clear  as  possible,  we  will  treat  only 
the  case  when  the  planet’s  potential  is  of  the  form: 


9  January  1963 


2 


SP- 1067/000/00 


It  will  be  clear  to  the  reader  that  everything  said  will  also  apply  to  the 
more  general  potential: 

(2)  V  -  (|i/r)  [1  +pE2BpPp(ain  p)/rP] 

In  equations  (1)  and  (2),  I  have  used  the  notation  of  Brouwer  (1959), 
namely,  (X,  k2,  and  Bp  are  physical  constants.  P  is  the  latitude,  r  the 
distance  from  the  center  of  the  planet.  P^sin  p)  are  Legendre  polynomials. 

In  terms  of  the  Delaunay  variables,  the  Hamiltonian  corresponding  to  the 
potential  (1)  can  be  written  (see  Brouwer,  1959): 


(3) 

F  “  F0(xi)  +  k2  Fll(xl'X2*  +  k2  Fi2(xi,x2,yl,y2) 

with 

(4) 

Fo  -  (ty^Xjx/Xj)2 

(5) 

F11  "  ("l^/2)  (1"3(x3/x2)2)  Ox/x^)3 

(6) 

F12  “  C°8  Jyl  +  fa  Qj<Xl,X2,*3>  COi 

(Jy!  +  2y2) 

where 

1 

2 

*  Gxa)  a  -  semi -major  axis 

1 

2  2 

x2  *  x^(l-e  )  e  ■  eccentricity 

x^  *  x2  cos  I  I  -  inclination 

I  -  inclination 
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*  mean  anomaly 


y^  =  argument  of  the  pericenter 


y^  ■  longitude  of  the  ascending  node 


(x^x^x^y^y^y^  are  canonical  variables  in  the  Hamiltonian  (3)).  See 
Poincar^  §  8. 

x^  is  a  constant  in  the  above  formulas  since  the  conjugate  momenta  y^  is 
not  present.  After  have  been  found  in  the  above  Hamiltonian,  one 
can  solve  for  y^  by  the  formula: 


(7) 


t 

y3  “  y3  'fQ  °VWyl'y2>dt- 


Since  (7)  is  merely  an  integration,  we  will  restrict  ourselves  to  solving 
the  problem  (3)  with  two  degrees  of  freedom.  Thus,  when  we  speak  of  periodic 
solutions,  we  mean  with  respect  to  the  four  variables  x^,x^,y^,y2«  It  will  be 
clear  to  the  reader  that  our  methods  would  also  carry  over  with  the  addition  of 
x^  and  y^;  but, once  again, to  increase  the  clarity  of  presentation  this  added 
complication  will  not  be  introduced.  Now  define 


(8) 


nl  *  SF0<xi>/dxi 


Any  quantity  throughout  this  paper  with  a  superscript  o  will  mean  it  is 
a  constant;  thus  y°  is  a  constant. 
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■  -  dFn(x°)/dXj 

If  n*  =  0,  we  have  the  so-called  critical  angle  case.  If  n*  /  0,  we  have 
the  non-critical  angle  case. 

Existence  of  Periodic  Orbits  in  General 

In  this  section  we  will  briefly  review  the  results  of  Poincare  (1892)  §  42 
on  the  existence  of  periodic  solutions.  We  will  show  that  these  results  are  not 
directly  applicable  to  the  Hamiltonian  arising  from  the  motion  of  a  satellite 
of  an  oblate  planet.  In  latter  sections  we  will  show  how  the  Hamiltonian  can 
be  transformed  so  that  the  Poincare  criteria  are  applicable. 

Poincare  (1892)  §  42  considers  a  Hamiltonian  of  the  form; 

(9)  H  -  H0(aj,a2)  WV 

where  a]*a2'vi'v2  are  canonical  variables,  corresponding  to  our  previous  x's 
and  yfs  with 

(10)  Ht  *  Z  B(a1,a2)  cos  (m^  +  4-  h(a]/a2)) 

and  where  the  summation  is  over  all  integers  m^,m2. 

If  we  set  €=0  in  (9)  and  solve  the  resulting  equations  of  motion,  we  obtain 
the  unperturbed  solution: 


a 


i 


(11a) 


o 

ai 
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(lib)  wi  =  "i  + 

(11c)  "j  =  "  ^Ho*ai^aj 

Poincare  further  assumes  the  a?  have  been  so  chosen  that  (n^/n^)  =  (p/q) 
with  p  and  q  integers.  Then  considering  w^,W£  angle-like  variables,  Poincare 
calls  the  unperturbed  solution  periodic  of  period  T  =  2irq/n^  *  2*p/n^.  p  may  be 
zero,  then  T  =  2rt/n^. 

The  principal  problem  is  to  find  periodic  solutions  of  the  same  period  T,  for 
small  values  of  €,  that  will  approach  the  unperturbed  solution  as  €  *>  o. 
Mathematically  this  can  be  formulated  as  follows.  Let 

(12a)  ai(0)  -  a°  +  bjL(e);  w.(0)  =  d°  +  c^e),  and 

(12b)  a.(T)  =  a°  +  b^e)  +  w.(T)  =  d°  +  c.(e)  +  nj  T  +  ^  . 

To  Insure  periodicity  and  the  desired  behavior  as  €  -►  o,  one  must  determine 
the  a°,d°,  80  that  ^  ■  0^  ■  o  for  all  sufficiently  small  €,  and  such  that 
b^(c),  c^(c)  approach  zero  with  €.  (These  are  exactly  the  same  equations  that 
Lefschetz  (1957)  p.  163  considers.) 

We  first  point  out  that  since  H  =  constant  is  an  integral  with  (dH/da^)  f  0, 
it  follows  that  if 

(13)  ^ ^2^ cl* *  ^1^1,^2,C1,€^  m  ^2^1,^>2*,C1^€^  ™  ® 

then  ^  will  also  be  zero  (see  Siegel  (1956)  top  of  page  126  or  Poincare, 
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Vol.  I,  p.  87).  Hence,  we  need  solve  only  (13).  Moreover  by  a  change  in  t, 
may  assume  =  c^(e)  ■  0. 

For  <f>^  one  finds 


(14) 


\-f  ((dWj/dO-n^dt  -  Ja-b H/da^-n^dt 
-  T((dHo(aJ+bi)/da1)-nJ)  +  0(e) 


similarly  for  ^ 


(15) 


♦2  “  T((dHo(a°+b1)/da2)-nJ)  +  0(e) 


For  ^  one  finds 


(16) 


i  i 

"  /«ydt)dt  »  f  (dH/dw2) 

T 

=  p/  (dH./dw2)dt  -  e^fHjJ/Sd®)  +  0(e2) 


where  if  is  defined  in  (10)  then 


(17)  [Hj]  -  T  Z’B(a°,a°)  cos  (n^d®  +  m2d°  +  h(aj,a°)) 


and  where  £' means  only  sum  over  the  terms  where 


(17a) 


mlnl  +  m2H2  "  0 


(This  follows  from  (16)  on  substituting  the  unperturbed  solution  in  H^.  Note 
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also  that  we  assume  d°  *  0. ) 

It  now  follows  from  the  implicit  function  theorem,  that  if  one  is  to 
solve  the  three  equations  (13)  for  small  values  of  €,  the  following 
conditions  are  sufficient: 


/  d2H 

d2H  /  d2H  \  2\ 

(18) 

f  o 

P  (  .  o.  )  ) 

da2da2  \da^da2J  J 

0 

ai =  ai 

and  d° 

*  is  so  chosen 

that: 

(19a) 

3fHl(a°,a°,^)]/ad°  =  0 

(19b) 

d2fH1(a®,a°,d^)  l/dd^dd®  +  0 

If  the  Poincare  criteria  (18)  and  (19)  are  fulfilled,  we  are  guaranteed 
periodic  solutions  for  sufficiently  small  values  of  e.  If  [H^]  ^  constant, 
since  [H^]  is  periodic  it  will  always  have  a  maximum  and  minimum;  thus  there 
will  always  exist  a  d^,  satisfying  (19a),  (19b). 

Let  us  now  see  where  the  Hamiltonian  (3)  fails  to  fit  these  criteria.  The 
obvious  choice  for  Hq,H^  is  (using  the  definitions  (3),  (4),  (5),  (8))  (we 
clearly  equate  with  a^,  and  y^  with  w^): 


(20) 


F  + 


k  F 
*2  11' 


eH, 


eF 


12 


where  in  (20)  we  consider  as  a  fixed  quantity,  €  as  variable.  Thus  (see  (11c)) 
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With  this  definition  of  Hq,  it  follows  that  the  development  in  the  present 
section  is  applicable.  Furthermore,  it  follows  that  for  most  values  of 
a°,a°,a°  (18)  is  fulfilled;  for  those  values  where  (18)  is  not  fulfilled  one 
can  treat  the  square  of  the  Hamiltonian  as  suggested  by  Poincare  §  43  to  obtain 
a  new  Hamiltonian  where  (18)  is  fulfilled.  Hence  (18)  causes  no  trouble. 

It  is  (19b)  that  is  not  fulfilled.  From  (6)  it  follows  that  m2  is  either 
0  or  2;  from  (21)  it  follows  that  n^  is  a  very  small  number.  Hence  there  are 
no  terms  in  H^,  where  (17a)  is  fulfilled;  thus  [H^]  =  0,  and  (19b)  is  not 
fulfilled. 

In  the  next  section  we  will  introduce  a  transformation  that  will  convert 
the  Hamiltonian  (3)  to  one  where  the  criteria  (18),  (19)  are  fulfilled. 

The  Poincare — von-Zeipel  Transformation 

It  is  clear  that  if  [F^]  =  0,  perhaps  if  we  develop  in  powers  of  €  we 
will  eventually  come  to  a  term  corresponding  to  [F^]  that  will  not  be  Identically 
zero.  The  most  straightforward  way  to  do  this  is  by  what  is  now  known  as  the 
von-Zeipel  transformation.  However,  it  should  be  mentioned  that  practically 
the  whole  of  Poincare  (1893)  Volume  II  is  devoted  to  this  same  type  of  trans¬ 
formation.  We  follow  principally  Poincare  §134  and  §125  in  our  development. 

Let 

ZP  "j-oeJsj('l'*2'yl'y2) 


(22) 
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with 

(23)  s0  =  a1y1  +  v2 

Now  consider  the  partial  differential  equation: 

(24)  F(3Ep/3y1^Ep/ay2,y1,y2)  =  C^a^  +  eC^a^a^  ... 

+  ePC?(ava2) 

where  F  is  the  Hamiltonian  (3)  (with  changed  to  e),  and  the  C^^Ca^^a^)  are 
to  be  determined. 

If  like  powers  of  €  are  equated  in  (24),  one  finds  the  series  of  equations 


(25)  P0Cal>  "  W 

(26)  n^SSj/ayj)  =  ^(SSj  j/Syj)  +  ^(SSj.j/Syj)  +  »j  -  C^.y.j) 

with  a  polynomial  of  the  j-th  degree  in  (dS^/dy^)  k  <  j-1,  the  coefficients  of 

this  polynomial  being  either  d^^a^/Sx"  or  dpF1(a^a2,y1,y2)/d<1Xjdp  **x2, 

and 


(27) 


To  any  solution  S 


(1) 

J 


o 


of  (26), 


aFo(a1)/Sx1 

^Fll(ai)/^x2 

aP12<ai,yl)/&*2 


we  can  always  add  on  a  solution  S 


(2) 

j 


that 


only  depends  on  y2« 


-  ■ 
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Let  us  first  proceed  formally  and  assume  at  some  stage  we  have  computed 

*  *  *Sj  2>  and  Sj^‘  We  n0W  USe  ecluation  (26)  to  comPute  and  Sj^# 

Since  is  assumed  known  we  have: 

q2(SS1-l/^y2)  +  ®i  =  5  A^al,a2^  COS  (mlyl  +  m2y2  +  h(al^a2)) 

■lf° 


(28) 


+  ZA^a^a^  cos  (m^j  +  h(aj,a2)) 

m2/° 


+  C'Ca^ap 
(2) 

Now  determine  C^a^a^)  and  to  eliminate  the  last  two  terms  on  the 

right-hand  side;  thus 


(29) 


°j<*l'*2>  "  C'<*1'*2> 


(30) 


S$2)  *  -2  (A,(a1,a2)/m2n2)  sin  (m2y2  +  h(a.,a2)) 
J  m.fo 


If  +  sj^2  is  used  In  (28)  rather  than  S^,  then  since  q* 

has  only  terms  which  depend  on  y^,  only  the  first  series  on  the  right-hand 
side  is  affected.  Thus  with  (28),  (29),  (30)  inserted  in  (26)  we  may  choose 
a  solution  of  the  form: 


(31) 


=  -  2  (A(a1,a2)/m1n°)  sin  (m^  +  m2y2  +  h(a1,a  )) 
m1r° 


With  this  method  of  iteration  we  determine  S^,...,S  .  Away  from  the 


o  1 


critical  angle  both  n^n  are  0(1);  hence  there  are  no  small  divisors  in  (30),  (31). 
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The  above  formal  method  of  proceeding  can  be  fully  justified.  First  we 
note  that  the  sum  in  (30)  contains  only  a  finite  number  of  terms;  hence  term 
by  term  integration  is  permissible. 

The  justification  of  term-by-term  integration  in  (31)  is  slightly  more 

3  3 

involved.  First  we  recall  that  essentially  =  (a/r)  +  (a/r)  cos  2(f+y2), 

and  then  the  following  formulas  (see  Brouwer  1959,  Wintner  1947): 

(a/r)  =  1+2  n|1Jn(ne)  cos  nyx 

2  00 

cos  f  *  -e  +  (1-e  )  E'J  .(ne)  cos  ny. 

n=  n  - 1  I 

sin  f  -  (1-e2)  1^_£|)Jn_1(ne)  sin  ny^^ 
dt/c^  =  (1/exp  (x2/x1)2(S\|f/3e) 
df/dx2  =  -(l/ex1)(x2/x1)(dt/Se) 
d(a/r)/de  =  (a/r)2  cos  f 
df/de  =  ((a/r)  +  (xj/x^2)  sin  f 

From  these  formulas,  together  with  the  estimates  for  J  (mz)  given  in 

m 

Wintner  1947  §294,  plus  standard  theorems  on  the  multiplication  and  integration 
of  Fourier  series,  it  follows  that  our  formal  procedure  is  justified. 

With  Ep,  Cj(a^,a^),  defined  as  above,  we  follow  Poincare 

in  changing  to  new  canonical  variables  a^,w^  by  the  formulas: 
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(32) 


xi  "  <asp/^i> 
wi  "  (Ssp/aai) 


The  Hamiltonian  becomes: 

(33)  H  =  F(x1#y1)  =  F(dEp/dyi)yi)  -  .  .+€°Cp+€P+1®t>fl(a1,a2,y1,y2,e) 

where  clearly  ®  is  periodic  of  period  2rt  with  respect  to  y^y^  anc*  expandable 
as  a  power  series  in  e. 

Write  the  change  of  variable  equations  as: 

(34)  xt  =  ai  +  S(£p-So)/3yi  *  wi  "  d(2p‘s0)/^*i 

Because  the  are  periodic  with  respect  to  y^y^,  follow  from  the 

above  equations  that  if  y ^  is  changed  to  y^  +  2^*  and  to  w^  +  2k^n}  the 

equations  will  not  change.  Hence  x^  and  y^-w^  are  periodic  of  period  2K  with 

respect  to  The  Hamiltonian  (33)  therefore  is  also  periodic  of  period 

2jt  with  respect  to  w^  and  w^.  Thus  the  Poincare--von-Zeipel  transformation  has 

the  very  important  and  interesting  property  that  it  takes  periodic  solutions  in 

a^w^  into  periodic  solutions  in  x^y^^  and  vice  versa.  Hence,  to  study  periodic 

solutions  it  does  not  matter  which  set  of  variables  we  use. 

In  summary:  Xj  and  y  are  functions  of  a^,w^,  being  periodic  of  period 

2n  with  respect  to  and  w^*  The  equations  of  motion  with  respect  to  the 

canonical  variables  a.,w.  comes  from  the  Hamiltonian: 

1'  1 


9  January  1963 


13 


SP- 106 7/000/00 


(35)  H  *=  HQ(a1,a2)  +  i.f+l^VW'VV 

P  i 

with  Ho(a1,a2)  -  F^)  +  ^(a^ap,  and  H.  of  the  form  (10). 

Periodic  Orbits  for  Non-critical  Angles 

We  note  that  we  started  with  an  of  the  form  (10),  where  had 

only  the  values  0  and  2.  Then  if  we  defined  [F^]  =  [H^]  as  in  (17),  [F^]  s  q. 
However,  in  the  construction  of  the  Poincare--von-Zeipel  transformation,  the 
series  (dS^/dy^)  are  multiplied  with  each  other  and  in  the  process  the  range  of 
values  of  grows.  Thus,  we  may  assume  that  for  non-critical  angles,  we  have 
transformed  to  (35)  and  that  we  have  chosen  p  so  large  that  [H^^]  ^  0* 
fact,  it  is  clear  from  our  method  of  construction  that  we  may  assume  [H  « 
A(al,a2)  008  ((P+l)^  +  ^(ax^a2^'  Hence,  the  Poincare  theory  of  periodic 
orbits  is  now  applicable  to  non-critical  angles  orbits,  with  n*j  defined  by  (35) 
and  (11c). 

Hence,  we  conclude  that  if  we  assume  the  initial  eccentricity  is  not  zero, 
(so  that  the  Pj,Qj  (®)  Hoes  not  vanish)  and  the  initial  inclination  angle  is 
not  the  critical  angle,  then  for  any  initial  value  of  semi-major  axis, 
eccentricity,  inclination  and  for  certain  values  of  the  initial  value  of 
perigee  (those  values  satisfying  the  Poincare  criteria  (19)),  the  equations  of 
motion  of  a  satellite  of  an  oblate  planet  will  have  periodic  solutions.  Note 
that  the  period  for  non-critical  inclinations  depends  on  the  period  it  takes 
the  perigee  to  make  one  revolution. 

Furthermore,  since  the  equation  giving  rise  to  the  initial  value  of 
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P+1 

perigee  comes  from  a  term  of  0(6^  )  where  p  is  very  large,  it  is  conjectured 
that  by  changing  the  oblate  planet 1  s  potential  very  slightly  it  would  be 
possible  to  assign  the  initial  value  of  perigee  arbitrarily  and  still  obtain 
a  periodic  solution. 

Periodic  Orbits  at  the  Critical  Angle 

At  the  critical  angle  n^  ■  0.  However,  irrespective  of  the  value  of  n^ 
if  we  apply  the  Poincare--von-Zeipel  transformation,  we  obtain  the  following 
equation  for  S^. 

n^PV^)  -  "  ci<ai>a2^ 

Because  contains  no  terms  of  the  form  cos  m2y2,  C^S  e<luat^on  *8  easily 
solved  (see  Brouwer  (1959)  p.  380-381).  Now  letting  Zp  ■  SQ  +  the 

transformed  Hamiltonian  is  of  the  form  (see  Brouwer  (1959)  p.  385): 

H  “  Val>  +  ePll*al'V  +  +  Al^al'*2^  COS  2w2 

q 

+  2  A2(a1>a2)  cos  (m^  +  ny^)]  +  0(e  ) 

m^o 

It  is  possible  to  apply  the  results  of  Poincare  §42  directly  to  this 
Hamiltonian.  The  unperturbed  Hamiltonian  is 

“o  ■  Fo<*l>  +  k2  P11<V*2> 


with 


V 
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«  -(aFo(a°)/Sa1)  -  k2  n(a°,a°) /ba^ 

*2  =  *k2  SFll(al'a2^a2  =  0 

thus  the  unperturbed  orbit  has  period  T  =  2it/n°  (there  is  nothing  in  Poincare 
§42  that  does  not  permit  n^  =  0;  in  fact,  in  §44  he  reduces  the  general  case 
to  this  situation). 

Now  [H  ]  AQ(a^,a2)  +  A^(a^,a2)  cos  2w2  and  thus  by  the  Poincare  criteria 

(18),  (19)  it  follows  that  at  the  critical  inclination  there  are  periodic 
orbits  of  period  T  *  2rt/n^  f°r  semi-major  axis  values,  and  all  eccentricities 
not  identically  zero,  and  for  initial  values  of  perigee  either  0  or  jt. 

(in  a  private  communication  W.  T.  Kyner  states  that  he  has  also 
obtained  this  result. ) 

Cone lusion 

The  existence  of  periodic  orbits  at  both  critical  and  non-critical  angles 
has  been  shown  for  all  semi-major  axis  and  eccentricity  values  and  for  the  positions 
of  perigee  fulfilling  (19).  The  frequency  of  the  periodic  orbit  at  the  critical 
angle  is  the  usual  anomalistic  frequency.  The  frequency  of  the  periodic  orbit 
at  non-critical  angles  is  the  frequency  of  the  motion  of  perigee  (draconic 
frequency) . 

An  interesting  problem  that  remains  to  be  discussed  is  how  the  periodic 
orbits  away  from  the  critical  angle  will  approach  the  periodic  orbits  at  the  critical 
angle  as  the  inclination  changes. 
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